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Introduction 

Computer-based  medical  modeling  and  simulation  offers  the  potential  for  a  completely 
new  way  to  learn  medical  care.  The  use  of  computers  as  human  assistants  for  medical  diagnosis 
and  training  will  require  new  knowledge  domains.  The  most  fundamental  among  these  is  the 
digital  characterization  of  biophysical  properties  of  perfused  organs  obtained  from  in  vivo 
measurements,  accounting  for  tissue  properties  integral  to  the  organ  as  well  as  boundary 
conditions  which  affect  the  organ’s  global  behavior. 

Our  research  program  is  designed  to  obtain  non-destructive  measurements  of  living 
tissue-liver,  spleen  and  kidney-from  large  animal  models  and  to  use  this  data  to  characterize 
non-linear  behavior  during  large  scale  deformations,  as  would  be  seen  in  clinical  practice.  We 
are  designing  instruments  to  measure  tissue  properties  and  then  creating  mathematical  models  to 
allow  representation  of  these  data  in  a  computer  simulation.  During  the  first  year  of  research,  we 
have  studied  small  scale  deformations  of  porcine  liver  in  vivo,  using  both  indentation  and  rotary 
shear  measurements.  We  developed  a  useful  standard  reference  object  to  examine  general 
mathematical  representations  of  deformation  and  published  this  standard  at  a  publicly  accessible 
website. 

All  work  has  been  done  according  to  institutional-  and  Defense  Department-approved 
animal  studies  protocols. 

Body 

The  following  sections  present  our  first  year  results  organized  according  to  the  Statement 
of  Work  from  our  original  research  proposal.  The  points  in  that  Statement  include: 

•  Design  whole  organ  measurement  tools:  indentation-shear  property  tool  and 
implanted  device 

•  Test  tools  and  obtain  measurements  on  ex  vivo  tissue 

•  Validate  measurements  using  hand-held  and  motorized  bench  tools  standards  for  ex 
vivo  tissue 

•  Initiate  FEM  models 

•  Begin  indentation  testing  of  whole  organs,  end  Year  1 

Design  whole  organ  measurement  tools 

Indentation-shear  property  tools 

Our  first  device  developed  to  measure  the  properties  of  tissue  is  the  Tissue  Material 
Property  Measurement  Tool  (TeMPeST  1-D),  which  employs  normal  (perpendicular)  indentation 
to  apply  forces  to  tissue  and  measure  the  force-displacement  response  [see  Appendix  A- 
Ottensmeyer,  2002a],  Since  the  ultimate  goal  of  our  research  is  the  measurement  of  human 
tissues,  this  early  instrument  was  designed  for  use  in  minimally  invasive  surgery.  Successful 
completion  of  the  instrument  has  demonstrated  that  building  such  an  instrument  is  possible,  and 
has  highlighted  numerous  areas  for  improvement  in  future  devices. 

The  TeMPeST  1-D  (see  Figure  1)  is  a  12mm  diameter  minimally  invasive  instrument, 
which  can  be  used  under  laparoscopic  conditions,  or  during  open  surgery.  It  can  generate 
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arbitrary  trajectories  with  a  range  of  motion  of  ±500pm,  forces  up  to  300mN,  and  bandwidth  of 
approximately  80Hz. 


Figure  1:  (left)  TeMPeST  1-D  instrument  passing  through  standard  1 2mm  surgical  port,  (right)  Cutaway  details  of 
assembled  sensor/actuator  package  and  details  of  moving  core  elements 

Through  a  collaboration  with  the  Minimally  Invasive  Surgery  Section  of  the  University 
of  Tubingen  (Germany),  we  were  also  able  to  make  use  of  the  Rotary  Shear  Applicator 
(ROSA2),  which  applies  a  torsional  shear  load  to  tissues,  while  measuring  the  angular 
displacement  of  the  contact  head.  The  device,  shown  in  Figure  2  (left  and  center),  fixes  the 
surface  of  the  tissue  with  an  outer  ring  (16mm  inside  diameter)  while  twisting  the  tissue  with  an 
inner  contact  head  (6mm  diameter).  The  issue  of  non-slip  contact  is  critical  for  proper  operation 
and  measurement,  so  a  variety  of  contact  heads  were  investigated,  including  a  densely  packed 
needle  array,  a  flat  tip  combined  with  surgical  adhesive,  and  a  novel  surface  modeled  after  the 
minute  “teeth”  of  sharkskin. 

The  ROSA2  does  not  exhibit  the  same  mechanical  response  as  the  TeMPeST  1-D,  so  a 
closed  loop  control  system  was  developed  for  it,  while  still  employing  the  same  interface 
hardware  and  software  as  the  TeMPeST  1-D. 

Implanted  device 

The  goal  of  the  implantable  device  is  to  measure  force,  displacement,  and  velocity  data 
on  arbitrary  local  volumes  within  a  whole  organ  to  determine  the  complex  material  properties  of 
the  organ’s  native  parenchymal  tissue.  A  proof-of-concept  device  to  accomplish  this  has  evolved 
from  the  percutaneous  retraction  device  described  in  the  initial  proposal.  The  idea  behind  the 
percutaneous  retraction  device  was  that  a  needle  could  be  inserted  into  an  organ  and  then  a  “T” 
would  be  deployed  that  would  move  within  the  organ  as  force,  displacement,  and  velocity  data 
are  collected.  However,  this  particular  device  is  better  suited  for  retracting  hollow  organs,  and 
does  not  allow  complete  control  of  device  configuration  within  the  tissue. 
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Figure  2:  (left)  R0SA2  in  contact  with  silicone  gel  test  object  (outer  ring,  optical  position  sensor  visible);  (center) 
R0SA2  with  outer  ring  removed,  showing  flat  tip;  (right)  detail  of  "rotary  sharkskin  "  contact  tip  design1 

The  T-needles  are  an  initial  prototype  designed  to  completely  control  the  device’s 
orientation  within  the  organ’s  tissue  (collinear  axis,  constant  initial  separation  between  the 
parallel  compression  bars,  perpendicular  “T”  configuration,  etc.)  such  that  a  uniform  strain  field 
is  applied  to  a  local  internal  volume  via  parallel  compression  bam  while  force  and  displacement 
data  are  measured  (Figure  3).  FEM  (finite  element  modeling)  aided  the  design  of  the  prototype 
by  determining  the  necessary  separation  between  the  bars  to  neglect  boundary  condition  effects 
(i.e.  area  under  investigation  is  small  compared  to  area  of  excised  tissue)  given  a  known  ex  vivo 
geometry  of  sample  tissue.  The  needles  are  aligned  using  the  v-block  set-up  and  then  two  side 
needles  containing  the  compression  bars  are  inserted  into  the  sample.  The  shafts  are  then 
inserted  and  screwed  into  the  compression  bars  before  the  ends  of  the  compression  needles  are 
removed  (see  Figure  4).  The  initial  device  is  manually  driven  using  a  micrometer  drive  plate 
while  force  data  is  collected  using  a  strain  gauge. 

Test  tools  and  obtain  measurements  on  ex-vivo  tissue 

Indentation-shear  property  tools 

A  variety  of  mechanical  elements  were  used  to  validate  the  performance  of  the  TeMPeST 
1-D,  including  mechanical  springs,  masses,  and  three  grades  of  silicone  gel.  Each  of  these 
materials  was  tested  independently  to  determine  their  properties,  so  that  the  values  measured  by 
the  TeMPeST  1-D  would  have  a  gold  standard  for  comparison.  The  details  of  these  tests  are 
described  in  [see  Appendix  A-Ottensmeyer,  2002b],  and  demonstrate  good  agreement  between 
TeMPeST  1-D  measurements  and  the  respective  standards. 

ROSA2  is  an  improved  version  of  an  earlier,  larger  device  which  was  subject  to 
excessive  friction  in  the  actuator.  The  initial  work  with  ROSA2  involved  determination  of  the 
calibration  curve  for  the  position  sensor2  and  the  torque  constant  for  the  galvanometer  motor. 
Simultaneously,  the  system  model  for  the  instrument  was  determined  so  that  a  digital  closed  loop 
controller  could  be  generated. 


1  This  device  has  a  diameter  of  6mm,  with  feature  sizes  on  the  order  of  0. 1  mm.  One  manufacturing  technique  (3  -D 
printing,  Z-Corp)  has  been  determined  to  be  unable  to  generate  this  geometry,  and  a  search  is  ongoing  for  other 
suitable  methods, 

2  An  optical  device  which  uses  two  polarized  filters  rotating  with  respect  to  each  other.  The  nature  of  this  system 
results  in  a  non-linear  (sinusoidal)  relationship  between  angular  position  and  sensor  output  voltage.  Calibration  and 
curve  fitting  permits  extraction  of  angular  data  from  the  raw  voltage. 
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T-Needle  compression  bar 


Figure  3:  T-needle  concept  and  initial  prototype .  (A)  T-needle  concept  showing  volume  of  parenchymal  tissue 
under  a  nearly  uniform  strain  field  from  displaced  T-needle  compression  bars.  The  strain  field  produced  goes  from 
highest  (red)  in  between  the  bars  to  lowest  (yellow)  outside  the  region  of  interest.  Force,  F(t),  and  displacement , 
x(t),  as  a  function  of  time  is  applied  and/or  measured  (B)  Top  view  and  (C)  angled  view  of  T-needle  prototype  set - 
up.  (1)  3mm  diameter  compression  bars  set  1.9  cm  apart,  (2)  3mm  diameter  shafts  that  move  compression  bars.  (3) 
Micrometer  drive  plate.  (4)  Strain  gauge  one  -axis  force  sensor.  (5)  V-block  support  and  alignment  system. 


Figure  4:  F reshly  excised  porcine  liver  test  with  T-needles.  (left)  Insertion  ofT-needles  into  vacuum  mounted  tissue 
sample,  (right)  Tissue  sample  ready  for  testing 

In  parallel  with  testing  of  the  instrument,  contact  heads  were  manufactured  and  tested  for 
non-slip  contact  with  various  materials  including  silicone  gels  and  ex  vivo  bovine  liver  samples. 
To  facilitate  this  testing,  the  apparatus  shown  in  Figure  5  was  created.  Rotation  can  be  held 
fixed  and  the  sample  observed  to  determine  at  what  angle  slippage  begins  to  occur.  For  both  the 
tissue  adhesive  (Johnson  and  Johnson  liquid  bandage)  and  the  needle  array,  angles  without  slip 
of  up  to  90°  were  achieved  under  static  conditions. 


t 


Figure  5:  (left)  static  test  rig  with  interchangeable  contact  heads;  (right)  test  box  with  bovine  liver  sample. 

Further  testing  was  done  on  bovine  liver  samples  to  verify  the  qualitative  results  of  finite 
element  analysis  (FEA)  done  at  Tubingen  with  regards  to  rotation  of  tissue  within  the  bulk  of  the 
organ.  To  achieve  this,  narrow  gauge  spinal  needles  were  used  to  inject  dilute  barium  sulfate 
solution  into  liver  tissue  and  the  static  test  rig  was  x-rayed  with  a  medical  fluoroscope  as  the 
contact  head  was  rotated  (see  Figure  6).  Numerical  analysis  of  the  images  is  ongoing,  but  initial 
results  show  agreement  with  the  FEA. 


Figure  6:  (left)  static  test  rig  with  spinal  needle  for  deposition  of  barium  sulfate  solution;  (right)  representative 
fluoroscope  image  showing  injection  paths  and  contact  head. 

Implanted  device 

The  initial  T-needle  prototype  was  manipulated  by  hand  and  was  tested  on  ex  vivo 
samples  of  freshly  excised  and  butcher  bought  porcine  liver  as  well  as  silicone  rubber  (Figure  4). 
From  these  ex  vivo  results  we  were  able  to  gain  some  insight  on  the  behavior  of  soft  tissue 
parenchyma  (clearly  a  viscoelastic  response).  These  initial  experiments  revealed  changes  that 
need  to  be  made  to  obtain  more  accurate  data.  Such  modifications  and  optimizations  are 
discussed  in  the  conclusions  section. 

Validate  measurements 

Indentation-shear  property  tools 

The  utility  of  the  TeMPeST  1-D  had  been  demonstrated  prior  to  die  start  of  this  work,  so 
the  validation  work  emphasized  the  comparison  between  two  testing  paradigms  (i.e.  indentation 
vs.  torsional  shear)  (see  Figure  7). 
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Figure  7:  Two  modes  of  tissue  deformation:  (left)  normal  indentation  -with  applied  force  (F)  and  displacement  (x) 
and  (right)  torsional  shear  with  respect  to  fixed  outer  ring,  showing  applied  torque  (t)  and  angular  displacement 
(9).  Characteristic  geometry  defined  by  indenter  radius,  a. 

Direct  comparisons  between  the  TeMPeST  and  ROSA2  were  performed  on  silicone  gel 
samples  and  on  porcine  liver  in  vivo.  Analysis  of  the  gel  samples  show  that  measurements  of  gel 
elasticity  are  the  same  to  within  better  than  a  factor  of  two  (see  Figure  8).  An  early  comparison 
of  the  same  tests  on  bovine  liver  samples  is  also  shown,  which  demonstrates  identical  frequency 
dependence  in  the  compliance  and  phase  data3. 


Frequency  (Hz)  Frequency  (Hz) 

Figure  8:  (left)  silicone  gel  sample  elasticity  measurements ;  (right)  bovine  liver  sample  compliance  measurements 
(units:  m/N for  TeMPeST \  V fV ]  for  ROSA2,  where  V$  is  position  sensor  output  and  Vt  is  the  drive  current  control 
signal ,  both  in  Volts) 

Implanted  device 

Each  measurement  made  on  the  organ’s  components  needs  to  be  validated  against  some 
physical  standard.  With  the  material  parameters  and  constitutive  law  defined  and  compiled  for 
each  component,  the  whole  organ  should  also  be  validated  against  a  physical  standard.  Lastly, 
due  to  the  trade-offs  made  between  fast  computation  time  and  accurate  displacements  in  real¬ 
time  computer  simulations  of  soft  tissue,  there  is  a  need  for  physical  validation  of  these 
algorithms. 

Testing  the  various  components  in  multiple  configurations  and  using  inverse  FE  methods 
will  serve  as  our  means  of  validating  each  component  of  the  organ. 


3  The  compliance  curves  are  shifted  due  to  the  differences  in  output  unit  systems  from  the  two  instruments. 
Analysis  to  bring  the  data  to  a  common  format  is  ongoing. 
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Soft  Tissue  Reference  Standard 


As  a  first  step  toward  a  standard  tissue  atlas  for  particular  solid  organs,  we  designed  an 
experiment  to  establish  a  benchmark  for  mathematical  modeling  of  tissue  deformations  under 
non-linear  displacements.  This  tissue  phantom  data  has  been  posted  on  our  website 
[http://biorobotics.harvard.edu/truthcube/]  and  can  be  used  as  a  reference  by  other  ^investigators 
who  wish  to  test  their  models  against  a  known  dataset.  r 

In  Year  1,  we  designed  and  fabricated  a  geometrically  regular  solid  test  organ  comprised 
of  a  well-defined  material.  This  “Truth  Cube”  is  an  8  x  8  x  8  cm  silicone  cube  with  discretely 
conspicuous  internal  fiducial  markers  which  can  be  precisely  tracked  during  computed 
tomographic  (CT)  imaging.  The  Cube  was  imaged  in  CT  under  known  loading  and  boundary 
conditions  (uniaxial  compression  and  spherical  indentation).  The  images  obtained  during 
loading  conditions  allowed  us  to  measure  internal  displacements  of  the  fiducials  and  create  a 
Finite  Element  model  of  large  scale  deformations  with  known  boundary  conditions  in  the  test 
object  [see  below,  Appendix  A-Kerdok  &  Howe,  2002], 

Initiate  FEM  models 

Finite  Element  Models  (FEM)  may  be  useful  in  medical  simulation  to  represent  complex 
systems  such  as  whole  organs.  They  are  generally  agreed  to  provide  good  representations  of  the 
mechanics  of  these  complex  systems.  Such  representations  can  be  used  in  device  design  (see 
Section  E.1,2)  as  well  as  material  property  prediction. 

For  comparison  and  validation  purposes,  we  created  FE  models  of  the  regular  geometric 
solid  organ  phantom  called  the  Truth  Cube,  [see  Appendix  A-Kerdock  et  al.,  2002],  The  FE 
models  represented  the  three-dimensional  system  under  uniform  compression  and  spherical 
indentation.  We  were  interested  in  determining  what  the  simplest  model  could  predict,  and  if 
more  robust  and  elaborate  FE  models  could  predict  the  more  complex  loading  conditions  that 
would  occur  under  larger  strain  and  global  deformations.  In  short,  is  finite  element  modeling  an 
accurate  method  of  representing  medically  important  organ  behaviors  in  surgical  simulations? 

We  found  that  the  finite  element  models  can  indeed  predict  the  behavior  of  the  material 
for  small  strains,  but  that  the  model  breaks  down  for  strains  greater  than  22%.  Soft  tissues  have 
highly  nonlinear  material  properties  and  are  subjected  to  large  deformations  that  exceed  30% 
strain  during  surgical  manipulation.  Therefore  we  concluded  that  although  FE  modeling  is  a 
useful  and  powerful  tool  for  modeling  small  strain  situations  of  many  materials,  it  is  merely  an 
approximation  method  whose  inherent  limitations  prevent  it  from  being  a  viable  whole-organ 
validation  method.  Physical  measurements  are  currently  much  better  suited  for  whole-organ 
validation.  As  part  of  our  continuing  research,  we  will  investigate  novel  methods  of 
mathematically  modeling  large  scale  deformations  in  our  out-year  program. 

Begin  indentation  testing  of  whole  organs,  end  Year  1 

In  vivo  testing  of  TeMPeST  ID  and  ROSA  2  on  the  same  organs  in  the  same  location  has 
begun  (see  Figure  9)  and  will  continue  into  the  next  year.  Analysis  will  be  performed  as  data 
becomes  available.  In  addition  to  measurements  of  tissue  compliance,  additional  evaluation  of 
the  non-slip  behavior  of  the  contact  heads  has  been  performed,  as  has  histological  analysis  of  the 
tissue  underneath  the  contact  heads.  The  former  tests  involved  performing  high-velocity  step 
commands  to  the  angular  position  of  the  ROSA2,  and  examining  the  torque-displacement 
response.  In  a  number  of  cases  with  the  needle  array,  evidence  of  slipping  was  observed  (e.g. 
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changes  in  the  mean  angular  position  of  the  head,  changes  in  angular  position  without 
corresponding  changes  in  applied  torque.)  No  slipping  was  observed  when  the  tissue  adhesive 
was  used  to  ensure  contact.  Histological  observations  showed  no  damage  occurred  during 
testing.  Angular  displacements  with  the  ROSA2  instrument  are  limited  to  ±15°,  which  is 
significantly  smaller  than  the  range  tested  statically,  so  eventual  testing  in  humans  should 
similarly  not  result  in  tissue  damage. 


Figure  9:  ROSA2  instrument  being  used  to  test  porcine  liver . 

Key  Research  Accomplishments 

•  Established  collaboration  with  U.  Tubingen  research  group  to  acquire  access  to  additional 
test  instrument 

•  Validated  comparability  of  torsional  shear  with  normal  indentation  as  a  reliable  method 
for  measuring  tissue  properties 

•  Initiated  parallel  tests  of  ROSA2  and  TeMPeST  1-D  instruments  on  porcine  liver  in  vivo 

•  Designed,  developed,  and  initially  tested  a  proof-of-concept  prototype  to  internally 
deform  a  local  volume  of  parenchymal  tissue  while  recording  force  and  displacement 
data 

•  Established  a  technique  to  create  a  physical  standard  of  3D  whole-organ  internal 
displacement  fields 

•  Developed  a  web  site  and  posted  the  initial  raw  data  of  the  3D  internal  strain  fields  of  the 
silicone  organ  phantom  (the  Truth  Cube) 

•  Showed  that  FE  Models  are  an  adequate  validation  tool  for  linearly  elastic  materials 
under  low  strains.  A  physical  validation  is  needed  for  strains  greater  than  approximately 
25%. 


Reportable  Outcomes 

Journal  Papers 

•  Kerdok  AE,  Cotin  SM,  Ottensmeyer  MP,  Galea  AM,  Howe  RD,  Dawson  SL:  Truth 
cube:  Establishing  physical  standards  for  real  time  soft  tissue  simulation.  Medical  Image 
Analysis  (submitted  for  publication). 
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•  Ottensmeyer  MP:  TeMPeST  1-D:  an  instrument  for  measuring  solid  organ  soft  tissue 
properties,  Experimental  Techniques  2002;  26,  48-50. 

Conference  Papers 

•  Ottensmeyer  MP,  Salisbury,  Jr.  JK,:  In  Vivo  Data  Acquisition  Instrument  For  Solid 
Organ  Mechanical  Property  Measurement,  MICCAI  2001.  Proceedings  of  the  Medical 
Image  Computing  and  Computer-Assisted  Intervention  4th  International  Conference, 
Utrecht,  The  Netherlands,  pp975-82.  14-17  Oct  2001. 

•  Ottensmeyer  MP:  In  vivo  measurement  of  solid  organ  mechanical  tissue  properties. 
Society  for  Experimental  Mechanics  Annual  Meeting,  Milwaukee,  WI.  pp328-333.  10- 
12  Jun  2002. 

•  Ottensmeyer  MP:  In  vivo  measurement  of  solid  organ  visco-elastic  properties. 
MMVR02/10.  Proceedings  of  Medicine  Meets  Virtual  Reality  02/10,  Newport  Beach, 
CA.  IOS  Press.  pp328-33.  23-26  Jan  2002. 

Invited  LecturesAA/orkshops/Meetings 

•  Ottensmeyer  MP,  Kalanovic  D,  Gross  J:  Comparison  of  indentation  and  rotary  shear  as 
modes  for  interrogating  soft  tissue  visco-elasticity,  SMIT2002,  Proceedings  of  the  14th 
Annual  Scientific  Meeting,  Oslo,  Norway,  5-7  Sept  2002. 

•  Kerdok,  AE,  Howe,  RD:  Characterizing  Soft  Tissues  for  Surgical  Simulation:  Probing 
Parenchymal  Properties,  MIT  Health  Sciences  and  Technology  2002  Forum  Abstract 

•  In  vivo  measurement  techniques  for  determining  mechanical  responses  of  soft  tissues. 
Tufts  University  Mechanical  Engineering  Seminar  Series,  10  Apr  2002. 

•  Kerdok  AE,  Cotin  SM,  Ottensmeyer  MP,  Galea  AM,  Howe  RD,  Dawson  SL:  Truth 
cube:  Establishing  physical  standards  for  real  time  soft  tissue  simulation,  International 
Workshop  on  Deformable  Modeling  and  Soft  Tissue  Simulation,  Bonn,  Germany,  14-15 
Nov  2001. 

Online  resources 

•  Truth  Cube  website:  http://biorobotics.harvard.edu/truthcube/ 

Conclusions 

Analysis  of  the  in  vivo  testing  is  ongoing,  but  qualitative  results  indicate  the  following: 
rotary  shear  does  not  cause  damage  to  the  tissue;  rotary  shear  as  a  deformation  mode  is  less 
sensitive  than  indentation  to  pulmonary  action;  maintaining  contact  with  tissue  using  adhesive  is 
more  reliable  than  the  use  of  a  needle  array  (which  may  become  clogged  and  slip);  but 
application  of  the  adhesive  is  an  additional  step  that  makes  the  measurement  protocol  more 
complex.  Finally,  a  closed-loop  control  algorithm  was  developed  to  support  use  of  the  ROSA2 
(which  is  not  as  amenable  to  open-loop  control  as  the  TeMPeST  1-D).  This  experience  will  be 
applied  to  the  TeMPeST  1-D  over  the  next  year  as  in  a  number  of  cases,  more  precise  control 
over  the  trajectory  generated  would  be  desirable.  In  addition,  future  generations  of  the 
TeMPeST  instrument  will  incorporate  design  changes  to  make  them  less  sensitive  to  pulmonary 
motions,  most  likely  by  grounding  themselves  against  (and  therefore  moving  with)  the  tissue 
rather  than  the  operating  table. 
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The  T-needles  are  an  initial  proof-of-concept  for  measuring  parenchymal  properties  of  ex 
vivo  tissue  samples.  However,  it  is  obvious  that  the  current  design  will  not  work  well  as  our  final 
“implantable”  in  vivo  device.  We  are  in  the  process  of  optimizing  the  T-needle’s  shape,  size  and 
performance.  The  initial  studies  revealed  a  need  to  reduce  the  stress  concentrations  and  friction 
in  the  needles  to  ensure  that  we  are  uniformly  compressing  the  tissue.  We  also  noted  the  need  to 
accurately  measure  the  viscoelastic  (time  dependant)  response  of  the  tissue,  and  so  the  device  is 
currently  being  equipped  with  a  motor  to  enable  the  application  of  fast  step  displacements  and 
loads  while  varying  strain  rate.  Lastly,  it  was  evident  that  in  vivo  or  whole  organs  in  a  controlled 
environment  would  give  much  more  repeatable  and  accurate  results  and  a  chamber  is  being  built 
that  will  perfuse  the  organ,  control  the  temperature,  and  minimize  blood  loss  during  testing. 

Once  we  have  accurate  and  reproducible  data,  we  will  work  to  redesign  the  system  to 
incorporate  multiple  force-displacement  modalities  (shear,  tension,  compression,  etc.)  and  to  be 
more  compatible  to  the  constraints  and  conditions  of  an  in  vivo  environment. 

Our  technique  for  recording  the  3D  internal  strain  fields  of  a  solid  volume  has  provided  a 
physical  standard  for  validation  of  whole-organ  measurements.  We  will  continue  this  work  by 
applying  the  technique  to  whole  organs  and  assessing  other  imaging  possibilities. 

FE  Modeling  will  continue  to  aid  in  our  device  design.  Once  more  organ  data  is 
collected  an  FE  model  of  each  test  configuration  will  be  developed  for  component  validation  and 
constitutive  law  determination. 

“So  What?” 

The  knowledge  gained  in  the  first  year  of  this  program  establishes  a  path  toward  in  vivo 
characterization  of  solid  organ  tissue  behavior  during  clinically  relevant  manipulations.  This 
information  will  be  crucial  to  creation  of  effective  medical  training  systems  which  can  be 
designed  to  reproduce  biological  properties  without  requiring,  for  example,  the  use  of  animals 
for  training.  Although  medical  simulation  holds  much  promise  for  changing  the  way  medical 
knowledge  has  been  imparted  for  millennia,  truly  effective  learning  of  what  to  expect,  how  to 
react,  and  how  to  correct  errors  is  the  essence  of  medical  wisdom.  Without  biological  properties, 
computer  graphics  are  just  cartoons,  devoid  of  relevance  to  patient  care.  But  with  biologically 
accurate  behavior,  computer  simulations  can  become  a  replacement  for  learning  on  humans  and 
animals.  This  work  establishes  one  pathway  to  obtain  the  biological  data  to  enable  this  transition. 
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Ottensmeyer,  2002a 

Ottensmeyer  MP:  TeMPeST  1-D:  an  instrument  for  measuring  solid  organ  soft  tissue 
properties,  Experimental  Techniques  2002;  26:  48-50. 

Ottensmeyer,  2002b  '■" 

Ottensmeyer  MP:  In  vivo  measurement  of  solid  organ  viseo-elastic  properties. 
MMVR02/10.  Proceedings  of  Medicine  Meets  Virtual  Reality  02/10,  Newport  Beach, 
CA.  IOS  Press.  pp328-33.  23-26  Jan  2002. 

Kerdok  &  Howe,  2002 

Kerdok,  AE,  Howe,  RD:  Characterizing  Soft  Tissues  for  Surgical  Simulation:  Probing 
Parenchymal  Properties,  MIT  Health  Sciences  and  Technology  2002  Forum  Abstract 

Kerdock  et  ai.,  2002 

Kerdok  AE,  Cotin  SM,  Ottensmeyer  MP,  Galea  AM,  Howe  RD,  Dawson  SL:  Truth 
cube:  Establishing  physical  standards  for  real  time  soft  tissue  simulation.  International 
Workshop  on  Deformable  Modeling  and  Soft  Tissue  Simulation,  Bonn,  Germany,  14-15 
Nov  2001. 
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FEATURE 


Biomechanics  Series:  Part  2 
by  M.P.  Ottensmeyer 


TeMPeST  l-D:  AN  INSTRUMENT  FOR  MEASURING 
SOLID  ORGAN  SOFT  TISSUE  PROPERTIES 


Surgical  simulation  is  a  field  that  shows  much  promise 
in  terms  of  permitting  surgeons  to  learn  new  skills, 
as  well  as  enabling  the  development  of  new  proce¬ 
dures  and  instruments  while  minimizing  the  use  of 
animal  subjects.  Essential  to  simulation  with  realistic  force 
feedback  are  material  property  data  for  the  wide  variety  of 
soft  tissues,  both  in  healthy,  and  diseased  conditions.  Prop¬ 
erties  are  known  to  change  significantly  after  death,  so  mea¬ 
surements  should  ideally  be  made  on  living  tissue.  Further, 
tissue  often  has  visco-elastic  behaviors,  so  measuring  prop¬ 
erties  over  a  range  of  frequencies  is  important.  This  paper 
presents  the  design  of  the  TeMPeST  l-D  (1-axis  Tissue  Ma¬ 
terial  Property  Sampling  Tool)  and  its  use  in  ongoing  mea¬ 
surements  of  the  visco-elastic  properties  of  solid  organ  tis¬ 
sues. 

As  with  engineering  materials,  the  mechanical  behavior  of 
living  tissue  is  governed  by  constitutive  relationships  that 
relate  stress  and  strain  (and 
often  strain  rate).  Part  of  the 
problem  of  modeling  tissue  is 
choosing  a  suitable  mathe¬ 
matical  expression  that  cap¬ 
tures  observed  behaviors  like 
elasticity,  stress  relaxation 
and  creep.  Once  a  model  is 
chosen,  then  the  particular 
parameters  that  characterize 
the  tissue  need  to  be  found. 

To  do  this,  one  needs  to  apply 
loads  or  displacements  to  tis¬ 
sue,  measure  the  responses 
and  then  take  geometry  into 
account  so  that  stress  and 
strain  (rate)  are  related  in¬ 
stead  of  force  and  displace¬ 
ment. 

DEFORMING  TISSUES 

Because  of  the  wide  variety  of  tissue  and  organ  types,  there 
are  many  different  ways  of  determining  tissue  responses. 
Porcine  brain  was  tested  by  Miller  et  al.1  using  a  spherically 
tipped  indenter  executing  a  position  ramp  and  hold  under 
computer  control.  Carter  et  al.2  have  developed  a  hand-held 
indenter,  which  has  been  used  on  human  liver  during  sur¬ 
gery.  This  instrument  measures  depth  of  indentation  by  com¬ 
paring  the  position  of  the  indenter  with  a  surrounding  ring 
that  rests  on  the  surface  of  the  tissue.  Tay  et  al.3  have  made 
use  of  a  small  three  degree  of  freedom  (DOF)  robot  arm 


Editor’s  Note:  The  new  Biomechanics  Feature  Series  will  provide  a  forum  for 
novel  experimental  techniques  in  the  areas  of  biological  and  biologically  inspired 
materials .  This  series  will  showcase  development  in  advanced  materials  systems 
with  an  emphasis  on  biological,  biomedical  and  biomimcry  applications .  Series 
Editor,  K.E.  Perry,  Echobio. 

M.P  Ottensmeyer  is  a  Lead  Investigator  for  Simulation  Group,  CIMIT,  Cambridge, 
MA 


equipped  with  a  six  DOF  force-torque  sensor  and  a  right 
cylindrical  tip  to  measure  porcine  esophagus  tissues  by  in¬ 
dentation  combined  with  shearing  (transverse)  motion.  A 
number  of  groups  have  modified  surgical  instruments,  equip¬ 
ping  them  with  force  and  position  sensors  (Scilingo  et  al.4), 
as  well  as  motors  (Brown  et  al.5)  to  support  computer  con¬ 
trol.  With  these  instruments,  the  response  of  tissue  to  grasp¬ 
ing  can  be  measured.  In  the  case  of  5,  the  grasper  jaws  are 
rectangular,  simplifying  the  conversion  from  the  force- 
displacement  curve  to  material  property.  At  least  two  groups 
have  developed  suction-based  techniques  in  which  a  tube  is 
sealed  against  tissue  and  the  pressure  lowered.  The  distance 
that  the  tissue  bulges  into  the  tube  is  measured  optically 
Aoki  et  al.6  use  a  clear  tube,  while  Kauer  et  al7  placed  an. 
angled  mirror  at  the  bottom  of  a  tube  with  an  offset  hole  cut 
in  the  end,  which  provides  a  profile  view  of  the  tissue  surface 
to  a  camera  placed  at  the  top  end  of  the  tube.  A  rotary  shear 
device  that  uses  a  circular  probe  with  a  dense  array  of  needle 

tips  to  rotate  the  tissue  while 
holding  the  surrounding  tis¬ 
sue  motionless  with  a  larger 
ring  (Kalanovic  and  Gross8) 
has  been  used  on  porcine 
liver. 

TeMPeST  l-D 
INSTRUMENT 

The  TeMPeST  l-D  (Ottens¬ 
meyer®)  (see  Fig.  1)  was  de¬ 
signed  explicitly  to  investi¬ 
gate  the  visco-elastic 
properties  of  tissue  under 
small  deformations.  As  test¬ 
ing  of  human  tissues  was  a 
long  term  goal,  a  minimally 
invasive  form  was  chosen  so 
that  it  could  be  used  either 
during  open  or  laparoscopic  surgical  settings.  Portability  was 
also  an  issue,  so  only  the  instrument  and  a  laptop  and  dock¬ 
ing  station  need  to  be  brought  into  the  operating  theatre. 
The  shaft  of  the  instrument,  which  contains  the  force  and 
position  sensors  and  a  voice  coil  actuator,  can  be  separated 
from  the  handle  so  that  the  shaft  can  be  sterilized.  When  in 
use,  the  instrument  is  rigidly  fixed  to  the  operating  table 
with  a  laparoscope  holder. 

The  instrument  has  a  range  of  motion  of  +/-500  pan,  can 
exert  forces  up  to  300  mN,  and  has  an  open  loop  mechanical 
bandwidth  of  approximately  80  Hz,  which  exceeds  the  band¬ 
width  of  human  hand-eye  and  reflex  control.  The  circular 
indenter  tip  has  a  diameter  of  5  mm,  and  the  overall  diam¬ 
eter  of  the  instrument  is  12  mm,  suitable  for  use  with  stan¬ 
dard  surgical  cannulas  (although  10  mm  and  5  mm  versions 
are  more  common).  Figure  2  shows  the  detail  of  the  sensor/ 
actuator  assembly,  including  the  compression  force  sensor 
(modified  Cooper  Instruments  &  Systems  (Warrenton,  VA) 


Fig.  I:  TeMPeST  l-D  instrument  inserted  through  12  mm 
surgical  cannula 
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Fig,  2:  Detail  of  sensor  /actuator  assembly.  Flexural  bearing- 
mounted  moving  core  (right)  consists  of  force  sensor  with 
circular  tip,  cylindrical  voice  coil  and  LVDT  position  sensor 
core.  Cutaway  view  (left)  shows  magnet  and  pole-piece 
assembly, 

LPM  562  500G),  the  custom  made  voice  coil  and  rare  earth 
magnet  assembly,  and  the  LVDT  position  sensor  (Schaevitz 
(Hampton,  VA)  099  XS-B). 

Performance  was  verified  on  a  variety  of  standard  objects 
and  materials  including  mechanical  springs,  standard 
masses  and  a  number  of  silicone  gels  (GE  RTV  silicone  gel 
6166).  Three  different  grades  of  the  latter  were  created  so 
that  the  results  from  a  parallel  plate  rheometer  could  be 
compared  with  the  TeMPeST  l-D  results. 

ANIMAL  TISSUE  TESTING 

The  instrument  has  been  used  to  measure  the  responses  of 
porcine  liver  and  spleen  in  vivo}  as  well  as  rat  liver,  kidney, 
lung,  heart  and  abdominal  muscle  in  vitro.f  For  the  in  vivo 
testing,  a  typical  animal  tissue  measurement  sequence  pro¬ 
ceeds  as  follows.  Ventilation  with  pure  oxygen  is  suspended 
briefly  so  that  breath  motions  are  not  measured  superim¬ 
posed  on  those  generated  by  the  instrument.  A  waveform, 
which  may  be  sinusoidal,  a  chirp  or  some  other  arbitrary 
signal  is  commanded  to  the  voice  coil  motor,  which  deforms 
the  tissue  with  the  5mm  punch.  Sampling  is  complete  in 
under  20  seconds,  at  which  time  ventilation  resumes.  A  typ¬ 
ical  view  through  the  laparoscope  during  testing  is  shown  in 
Figure  3,  Since  the  sampling  interval  is  short  and  pure  ox¬ 
ygen  is  used  for  ventilation,  no  harm  comes  to  the  animal 
during  testing.  Similarly,  since  the  displacements  of  the  tis¬ 
sue  surface  are  small,  no  damage  is  incurred  from  the  in¬ 
strument.  These  issues  are  essential  to  consider  early  on, 
especially  because  successor  instruments  are  planned  for  hu¬ 
man  use. 

TA11  animal  testing  was  approved  by  the  relevant  Institutional  Animal  Care  and 
Use  Committees. 


Fig.  3:  Laparoscopic  view  of  indenter  tip  in  contact  with 
porcine  liver 


While  the  results  of  testing  are  frequency  (and  load)  depen¬ 
dent  stiffness  values,  they  are  dependent  on  the  geometry  of 
the  indenter  and  the  tissue  being  tested.  For  large,  solid  or¬ 
gans,  a  useful  first  approximation  is  that  the  tissue  is  iso¬ 
tropic  and  homogeneous  (generally  valid  for  liver  or  spleen), 
incompressible  (approximately  true  for  water-rich  living  tis¬ 
sue)  and  linear.  This  last  assumption  depends  on  the  use  of 
small  deformations,  relative  to  the  size  of  the  organ.  Further, 
if  the  small  deformation  approximation  is  valid,  then  the 
organ  can  be  treated  as  a  semi-infinite  body,  which  provides 
for  closed  form  relationships  between  Young's  modulus  and 
the  measured  stiffness  and  the  geometry  of  the  indentation 
probe.  In  Equation  (1)  (adapted  from  Hayes  et  al.10),  k  is  the 
measured  stiffness,  a,  the  radius  of  the  indenter,  if,  a  geo¬ 
metric  factor  which  is  unity  for  a  semi-infinite  body  and 
greater  than  one  for  a  thin  material  on  a  rigid  substrate, 
and  finally  E  is  Young's  modulus  (all  in  any  consistent  set  of 
units). 


This  expression  can  be  extended  to  a  frequency  dependent 
form  if  k  and  E  are  replaced  by  functions  of  frequency  allow¬ 
ing  viscous  and  inertial  effects  to  be  taken  into  account.  An 
example  of  the  results,  showing  E  as  a  function  of  frequency 
and  nominal  stress  (preload  divided  by  indenter  area»  is 
shown  in  Figure  4. 

For  the  case  when  the  organ  is  small  compared  with  the 
motion  of  the  instrument  (e.g.  rat  kidney),  or  when  the  tissue 
model  is  more  complex  than  the  extremely  simplified  one 
described,  finite  element  models  can  be  employed,  iteratively 
approaching  the  observed  response  by  modifying  the  para¬ 
meters  given  the  measured  geometry  and  stiffness.  Analysis 
of  the  rat  kidney  and  liver  data  is  ongoing. 

SUMMARY 

Measuring  the  material  properties  of  living  tissue  is  a  com¬ 
plex  problem,  not  only  because  of  the  nature  of  the  tissues, 
but  also  because  of  the  limits  imposed  by  minimizing  or  pre¬ 
venting  damage  to  the  test  subject.  As  a  result  a  wide  variety 
of  techniques  have  been  and  continue  to  be  developed,  in- 
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Fig,  4:  Frequency  dependent  and  non-linear  elasticity  of 
porcine  liver.  Squares  are  results  from  sinusoidal  tests,  point 
clouds  from  chirp  response. 


eluding  non-invasive  ones,  employing  external  stimulation 
and  imaging  technology,  as  well  as  invasive  techniques 
which  directly  contact  and  deform  the  tissue  to  be  studied. 
No  one  technique  is  suitable  to  capture  the  full  range  of 
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behaviors,  including  non-linear  stiffness,  visco-elasticity, 
anisotropy  and  inhomogeneity  However  each  has  its  own 
advantages,  and  can  be  used  to  contribute  data  to  and  cross 
check  the  results  of  the  others. 
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To  support  the  ongoing  development  of  software-based  surgical 
simulation  systems,  work  is  underway  to  acquire  the  mechanical 
properties  of  living  tissue.  When  such  simulations  include  force 
feedback,  visco-elastic  properties  must  be  evaluated  over  a  range  of 
frequencies  relevant  to  human  perception  and  motor  control.  A 
minimally  invasive  instrument  has  been  developed  which  can  perform 
normal  indentation  on  solid  organs,  and  apply  and  measure 
deformations  over  a  frequency  range  from  DC  to  approximately 
100Hz.  Measurement  performance  was  validated  on  a  series  of 
objects  and  materials  with  known  properties,  and  the  device  was 
subsequently  used  in  in  vivo  tests  on  porcine  liver.  Results  of  these 
validation  tests  as  well  as  the  data  extracted  from  the  in  vivo 
experiments  are  presented.  Testing  is  ongoing,  and  will  be  expanded 
to  more  completely  characterize  liver,  as  well  as  porcine  spleen  and 
other  solid  organ  tissues.  While  these  animal  tissue  property  tests  are 
valuable  in  and  of  themselves,  they  pave  the  way  for  the  development 
of  instruments  and  experimental  protocols  suitable  for  the 
measurement  of  human  tissue  properties. 


1.  Introduction 

One  of  the  difficulties  in  developing  realistic  computer-based  simulators  for  surgery  is  the 
scarcity  of  material  properties  for  living  tissue.  Models  which  capture  the  mechanical 
behavior  of  tissues  typically  use  estimates  or  in  vitro  measurements  of  the  model 
parameters  for  given  types  of  healthy  or  diseased  tissues. 

While  estimates  from  experts  may  be  suitable  for  surgical  training,  new  instrument 
and  new  procedure  prototyping  will  require  more  precise  property  values  to  closely  match 
real  organs.  In  vitro  testing  is  useful  as  a  place  to  begin,  but  because  of  changes  like 
absence  of  blood  perfusion  and  oxygenation,  changes  in  boundary  conditions,  temperature, 
etc.,  one  must  also  pursue  in  vivo  acquisition  to  get  accurate  tissue  property  data. 

Researchers  are  pursuing  non-invasive  [1,  2],  minimally  invasive  [3,  4]  and  open 
surgical  techniques  [5,  6,  7].  Among  the  non-invasive  techniques,  are  ultrasound  [1]  and 
MRI  elastography  [2],  which  make  use  of  static  or  dynamic  deformations  at  the  surface  of 
the  body,  and  observations  of  the  strain  field  imposed  within  the  body.  By  performing 
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inverse  calculations,  elasticities  of  tissues  relative  to  each  other  can  be  determined.  This  is, 
however,  a  complex  process,  and  is  highly  dependent  on  knowledge  of  the  relevant 
boundary  conditions.  The  minimally  invasive  and  open  surgical  techniques  impose  local 
deformations  on  tissue,  and  examine  relatively  local  responses,  [3]  and  [4]  employ  surgical 
instruments  equipped  with  sensors  to  measure  the  local  force-displacement  response,  [5], 
[6]  and  [7]  use  larger  instruments  to  impose  tensile  loads,  indentation  or  suction  to  tissues. 
By  comparing  the  force-displacement  responses  with  either  simplified  or  finite  element 
models  of  the  tissue,  again  the  inverse  problem  is  solved  to  extract  the  tissue  properties. 
These  latter  techniques  can  be  computationally  less  expensive,  and  can  impose  larger 
deformations  on  tissue  than  the  non-invasive  techniques.  However,  they  suffer  limitations 
in  that  local  measurements  must  be  repeated  at  numerous  locations  over  a  given  organ  to 
yield  a  complete  set  of  data.  In  addition,  for  organs  with  capsules  (e.g.  liver,  spleen,  or 
kidney),  it  may  be  difficult  to  differentiate  the  contributions  to  the  measured  stiffness  of  the 
parenchyma  versus  the  capsule,  since  both  structures  are  deformed  simultaneously.  For 
liver  and  spleen,  which  have  relatively  thinner/weaker  capsules  than  kidney,  this  effect 
might  be  negligible,  but  further  experiments  will  be  necessary  to  find  out  what  those 
contributions  are. 

In  vivo  tissue  property  measurement  is  still  a  young  field,  and  the  broad  range  of 
tissue  types,  the  number  of  parameters  needed  to  describe  them,  and  the  issues  of  age,  type 
of  disease  and  individual  differences  mean  that  there  is  still  a  need  for  additional  data.  No 
single  technique  will  likely  provide  all  of  the  data  necessary,  so  all  reliable  data  is  a 
welcome  contribution. 


2.  Methods  and  instruments 

This  research  makes  use  of  a  single  axis  indentation  device  called  the  Tissue  Material 
Property  Sampling  Tool  (TeMPeST  1-D)  [8,  9,  10].  It  can  pass  through  a  12mm  surgical 
cannula  or  be  used  in  open  surgery,  and  generates  small  amplitude  oscillations  up  to 
approximately  100Hz,  Range  of  motion  is  +/-500um,  and  the  maximum  force  that  can  be 
applied  with  a  5mm  right  circular  punch  is  300mN,  Oscillations  are  generated  by  driving  a 
voice-coil  motor  with  a  voltage  controlled  current  amplifier. 


Figure  1.  (left)  Detail  of  sensor/actuator  tip  of  TeMPeST  1-D  and  moving  core  elements.  12mm 
overall  diameter.  Force  sensor  is  compression-only,  suitable  for  indentation  testing.  Flexural  bearings 
provide  friction  free  support  for  voice-coil  motor-driven  core  elements.  Non-contact  position  sensing 
performed  with  LVDT  sensor,  (center  &  right)  TeMPeST  1-D  in  validation  tests  on  springs  and  inertial  loads. 
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2.1.  Validation  Tests 


Figure  2,  Frequency  dependent  compliance  tests  on  springs  and  masses.  For  springs,  TeMPeST  1-D 
frequency  sweep  data  agree  well  with  0.1  Hz  fixed  frequency  (quasi-static)  data,  also  with  cartilage  press 
instrument.  For  masses,  fitting  inertial  load  characteristic  to  data,  yields  good  measurements  of  mass. 


To  verify  that  the  instrument  could  be  relied  upon  to  provide  accurate  data,  and  to 
determin  measurement  limitations,  a  series  of  tests  on  mechanical  springs  (compliant 
loads),  known  masses  (inertial  loads)  and  silicone  gels  (visco-elastic  materials)  were 
performed.  When  tested  over  a  range  of  frequencies,  springs  should  exhibit  a  constant 
compliance,  while  for  masses,  the  compliance  should  fall  by  two  orders  of  magnitude  while 
frequency  increases  by  one  order.  As  shown  in  Figure  2,  both  of  these  results  are  observed. 
Further,  for  the  springs,  the  TeMPeST  1-D  measurements  agree  with  those  made  using  a 
small  compression  testing  device  (cartilage  press).  Similarly,  the  masses  calculated  are 
close  to  the  masses  employed  during  testing. 

Springs  and  masses  have  simple  dynamic  characteristics,  so  a  series  of  silicone  gels 
exhibiting  visco-elastic  responses  within  the  bandwidth  of  the  TeMPeST  1-D  were  created 
and  tested.  The  gels  were  independently  tested  using  a  parallel  plate  rheometer,  and  again 
good  agreement  was  obtained  (see  Figure  3). 


Figure  3.  Visco-elastic  material  testing.  Three  different  mixtures  of  two-part  RTV  silicone  rubber 
were  made  and  tested  using  a  parallel  plate  rheometer  (symbols)  and  compared  with  chirp  data  (point  clouds) 
taken  with  the  TeMPeST  1-D.  For  softest  gel,  force  sensor  signal  to  noise  ratio  is  too  small  to  obtain  reliable 
data,  providing  lower  limit  of  tissue  stiffness  that  can  be  measured  without  modifications  to  the  instrument. 
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Figure  4,  TeMPeST  1-D  instrument  in  use  measuring  porcine  liver  in  vivo.  TeMPeST  1-D  visible  in 
foreground  (left),  mounted  on  flexible  laparoscope  holder  arm.  Testing  guided  under  laparoscopic  camera 
guidance  (handheld,  background,  and  view  shown  at  right). 


The  TeMPeST  1-D  is  being  used  in  ongoing  in  vivo  testing  on  porcine  solid  organs1 
to  investigate  viscous  and  non-linear  elastic  behaviors  over  the  range  relevant  for  haptics  in 
surgical  simulation.  Figure  4  is  an  image  from  early,  proof-of-coneept  testing,  confirming 
that  the  instrument  could  be  successfully  employed  in  a  minimally  invasive  environment. 
One  disadvantage  of  this  particular  prototype,  however,  is  that  given  the  convex  shape  of 
most  of  the  desirable  internal  organ  surfaces,  it  is  possible  to  make  normal  contact  at  at 
most  one  location  on  a  given  organ.  To  obtain  data  over  the  surface  of  the  organ,  one 
would  require  numerous  insertion  points,  at  least  when  used  minimally  invasively. 

A  series  of  tests  was  begun  in  conjunction  with  other  researchers  performing  open 
surgical  testing,  which  afforded  access  to  large  regions  of  the  porcine  liver.  Anesthesia  was 
achieved  using  an  initial  injection,  and  was  maintained  using  an  inhaled  anesthetic 
throughout  the  testing.  The  pig  was  placed  on  pure  oxygen  ventilation,  which  could  be 
suspended  for  periods  of  approximately  20  seconds.  During  this  time,  the  TeMPeST  1-D 
was  brought  into  contact  with  the  surface  of  the  organ  with  a  predetermined  pre-load  force; 
the  indentation  waveform  was  generated  with  die  voice-coil  motor  and  imposed  on  the 
tissue  for  approximately  16.5  seconds;  the  TeMPeST  1-D  was  removed  from  the  tissue  and 
ventilation  resumed.  Blood  oxygen  saturation  was  monitored  and  was  not  observed  to  fall 
below  95%  during  testing. 


3.  Results 

Tests  on  porcine  liver  have  been  performed  using  both  sinusoidal  and  chirp  signals, 
covering  the  frequency  range  of  0.01  to  100Hz,  and  preloads  from  8  to  75mN.  Within  this 
range,  the  elastic  modulus,  calculated  from  the  measured  stiffness  and  the  assumption  of  a 
semi-infinite,  incompressible  elastic  medium,  has  been  shown  to  be  a  non-linear  function  of 
mean  applied  stress  (and  therefore  with  strain)  and  frequency.  A  preliminary  fit  to  the 
fixed  frequency  sinusoidal  data  takes  the  form: 

l°gi0  i?  =  3.31 +0,161  logI0  /  +  2.15  E-  4a*med 

Equation  1 


1  All  in  vivo  testing  was  performed  under  protocols  approved  by  the  animal  use  committees  of  all  of 
the  institutions  involved.  To  minimize  the  number  of  animals  used,  testing  was  performed  either  as  an 
addendum  to  a  pre-existing  protocol,  or  was  performed  together  with  other  researchers  conducting  similar 
tests. 
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Figure  5,  Chirp  and  fixed  frequency  sinusoidal  data  obtained  in  vivo  from  porcine  liver. 

where  E  is  the  elastic  modulus  in  Pa, /is  frequency  in  Hz,  and  o*med  is  the  median 
applied  stress  (median  force/punch  area).  These  data,  and  the  fit  are  shown  in  Figure  5. 
What  is  clear  from  the  testing  performed  to  date  is  that  the  chirp  tests,  while  convenient  in 
terms  of  acquiring  a  range  of  frequency  dependent  data  in  one  trial,  have  poor  accuracy,  so 
further  testing  will  likely  be  performed  using  fixed  frequency  sinusoids. 

Recognizing  that  the  elastic  modulus  is  the  derivative  of  stress  with  respect  to  strain, 
this  expression  can  be  rewritten  as: 

ct  =  -2020  ln(l  -0.0005e7'62+ol61ta/e) 

Equation  2 

This  fit  is  based  only  on  the  form  exhibited  by  the  data,  rather  than  necessarily 
having  a  physical  basis.  As  more  data  become  available,  parameters  for  more  traditional 
models  will  be  extracted. 


4.  Conclusions  and  ongoing  work 

The  TeMPeST  1-D  instrument  for  measuring  visco-elastic  properties  of  solid  organ  tissues 
has  been  completed.  Validation  testing  on  known  materials  was  performed  to  verify  the 
device’s  performance,  and  it  has  been  used  to  gather  frequency  and  stress  dependent 
elasticity  data  for  porcine  liver. 

Testing  is  ongoing,  and  tests  planned  for  the  remainder  of  2001  and  early  2002 
include  investigations  of  variation  between  individuals,  inhomogeneity  across  organ 
surfaces,  and  tests  of  porcine  spleen. 

Based  on  the  learning  experience  of  these  teste,  the  testing  system  will  be  further 
developed  to  be  sensitive  to  a  wider  range  of  tissue  stiffnesses;  be  less  sensitive  to 
disturbances  caused  by  cardiac  action;  and  eventually  be  modified  for  use  on  human  organs. 
Ultimately,  such  properties  will  be  used  in  creating  surgical  simulation  systems,  which 
would  provide  realistic  force-feedback  to  end  users  from  surgical  students  learning  the 
basic  skill  required,  to  experienced  surgeons  refreshing  their  skills  or  developing  new 
procedures. 
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With  the  onset  of  minimally  invasive  surgery  (MIS)  surgeons  are  inquired  to  learn 
new  skills  for  manipulating  and  visualizing  3D  tissue  from  2D  images.  Real  time 
simulations  of  MIS  are  needed  for  training  and  planning .  Developing  realistic 
simulations  is  a  formidable  challenge .  The  overall  project  aim  is  to  improve  real 
time  surgical  simulation  by  developing  mathematical  models  of  soft  tissues  (liver, 
kidney ,  and  spleen)  under  slow  deformations  typical  of  surgical  manipulation , 
implementing  new  devices  for  measuring  the  necessary  soft  tissue  material 
parameters  for  these  models ,  and  determining  the  natural  boundary  conditions  of 
these  organs .  My  focus  is  to  develop  tools  needed  to  measure  the  material 
parameters  of  soft  tissue  parenchyma f  model  the  observed  behavior  in  the  form  of  a 
constitutive  law ;  and  validate  both  the  parenchyma  and  whole  organ  results. 

The  complex  nature  of  soft  tissue  constituents  (capsula,  parenchyma,  and  vessels) 
needs  to  be  addressed  to  make  meaningful  material  property  measurements.  Tissue 
testing  devices  must  account  for  tissue  inhomogeniety,  anisotropy,  geometric  and 
mechanical  nonlinearity  under  large  deformations,  incompressibility,  loading  rate 
dependency,  and  viscoelasticity.  A  device  has  been  developed  to  arbitrarily  probe 
an  internal,  local  volume  of  parenchyma  under  a  uniform  stress  distribution  ex  vivo. 
The  device  applies  known  displacements  while  measuring  forces  between  2 
cylindrical  bars.  Initial  studies  of  the  “T-Needle”  prototype  have  been  carried  out  on 
butcher  bovine  liver,  silicone  rubber,  and  freshly  excised  bovine  liver.  Results 
indicate  a  viscoelastic  response.  Modifications  to  the  device  are  needed  to  produce 
repeatable  results,  and  before  in  vivo  studies  can  be  performed. 

A  mathematical  model  will  be  developed  to  define  the  constitutive  law  of  the 
parenchyma.  Either  a  quasilinear  viscoelastic  type  model  with  a  strain  energy 
function  elastic  response  (i.e.  Mooney  Rivlin),  or  a  biphasic  model  will  be  used 
initially  to  model  the  tissue’s  force-displacement  response.  Estimated  material 
parameters,  test  loading  conditions,  and  the  measured  displacements  will  serve  as 
inputs  to  a  finite  element  model  that  will  output  a  predicted  force  to  compare  against 
the  measured  force.  The  chosen  law  will  be  optimized  by  testing  it  in  several 
different  loading  conditions  until  the  estimated  material  parameters  converge. 

The  optimized  model  for  each  soft  tissue  constituent  will  be  combined  to  model 
the  organ  of  interest.  A  new  technique  is  under  development  to  validate  this  overall 
model.  The  goal  is  to  provide  accurate  3D  internal  strain  field  measurements  that 
can  be  compared  to  the  model’s  results  under  known  boundary  conditions.  A  “test 
organ”  has  been  developed  using  a  silicone  rubber  model  with  internal  radiopaque 
markers  whose  positions  are  tracked  under  known  loading  conditions  using  CT.  A 
similar  approach  will  be  used  next  on  actual  organs  for  overall  validation. 

This  work  is  pan  of  i  project  in  conjunction  with  the  CIMIT  Simulation  group  at  Massachusetts  General 
Hospital,  funded  h\  '  s  Army  Medical  Research  and  Material  Command,  and  the  Whitaker  Foundation. 
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Abstract 

Accurate  and  real-time  models  of  soft  tissue  behavior 
are  key  elements  in  the  development  of  medical 
simulation  systems.  However,  precise  validation  of 
these  models  remains  a  challenge.  Currently,  real-time 
modeling  is  at  best  validated  against  FEM  models  that 
have  their  own  intrinsic  limitations.  This  study  is  the 
first  in  a  series  that  will  develop  physical  standards 
resulting  in  a  database  of  relevant  information  used  to 
validate  real-time  soft  tissue  deformation.  In  this  first 
study,  a  simple  8cm  cube  of  silicone  rubber  with  an 
exact  pattern  of  embedded  Teflon  beads  was  tested  in 
uniaxial  compression  while  CT  images  were  taken.  The 
known  material  properties,  geometry,  and  carefully 
controlled  boundary  conditions  resulted  in  a  complete 
set  of  volumetric  displacement  data.  The  results  were 
also  compared  to  an  FEM  analysis  of  an  identical 
situation.  This  work  has  served  as  a  proof  of  concept 
for  a  robust  physical  standard  for  use  in  validating  real¬ 
time  models.  A  web  site  has  been  created  to  provide 
access  to  our  database  at: 
http://biorobotics.harvard.edii/truthcube/, 

1.  Introduction 

Fast,  realistic  mechanical  modeling  of  soft  tissues 
remains  a  formidable  challenge  in  the  development  of 
new  medical  applications.  Medical  training,  surgical 
planning,  and  image-guidance  systems  require  real¬ 
time  calculation  of  tissue  shape  and  interaction  forces 
as  the  user  manipulates  the  simulated  tissue. 

Implementation  of  these  modeling  systems  has  proved 
difficult.  The  most  frequently  used  approaches  to 
simulating  tissue  deformation  in  real-time  have  been 
spring-mass  models  and,  more  recently,  finite  element 
methods  (FEM).  This  has  proved  relatively 
unsatisfactory:  surgical  applications  involve  large 
deformations  of  nonlinear  soft  materials  and  require 
very  rapid  calculations.  Spring-mass  systems  have 


limited  accuracy  but  work  in  real-time,  while  FEM  has 
proved  successful  for  modeling  small  deformations  of 
linear  engineering  materials  where  calculation  speed  is 
not  a  paramount  concern.  Many  schemes  have  been 
proposed  to  overcome  these  limitations  [e.g.  PicOl, 
DeOl,  Ast98,  Kuh97,  and  Cot96],  often  balancing 
execution  speed  and  mechanical  accuracy.  The  nature 
of  this  tradeoff  is  not  well  understood,  and  it  is  difficult 
to  evaluate  the  relative  performance  of  new  systems 
when  speed,  accuracy,  tissue  properties,  and  geometry 
all  vary. 

Progress  towards  effective  real-time  simulation  would 
be  aided  by  a  consistent  means  of  evaluating  tissue 
mechanics  models.  Ideally,  this  evaluation  method 
would  provide  a  “gold  standard”  of  experimental 
results  for  comparison  with  simulation  results.  This 
would  take  the  form  of  fully  specified  material 
properties  and  stress,  strain,  and  displacement  fields 
throughout  the  tissue,  under  a  range  of  fully  specified 
and  surgically  relevant  boundary  conditions.  The  output 
of  new  modeling  systems  could  then  be  compared  to 
these  standard  results  to  validate  and  benchmark  their 
performance. 

We  are  working  to  create  such  a  database  of 
experimental  measurements.  The  goal  is  to  incorporate 
measurements  for  rubber  phantoms  with  simple 
geometry  and  for  a  range  of  more  complex  materials 
including  biological  tissues  and  organs.  These  results 
will  be  archived  on  a  web  site  for  widespread  access. 
This  paper  reports  our  initial  effort,  on  the  deformation 
of  a  cube  of  soft  polymer.  This  8  cm  “truth  cube”  is 
made  of  a  well-characterized  silicone  rubber  with 
fiducials  embedded  in  a  1  cm  grid  pattern  throughout 
its  volume.  We  have  acquired  computer  tomography 
(CT)  scans  of  the  cube  in  undeformed  and  loaded 
states,  and  extracted  the  displacement  of  each  fiducial 
using  image  processing  techniques.  Here  we  report 
these  displacement  results,  along  with  details  of  the 
cube  construction  and  boundary  conditions  in  the 
loading  tests.  We  conclude  with  a  discussion  of  plans 
for  extension  of  this  work,  and  an  invitation  for 
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participation  by  others  interested  in  the  validation  of 
fast  large-deformation  tissue  modeling. 

2.  Methods 

We  have  selected  a  cube  as  a  simple,  regular  shape  for 
this  initial  effort  to  assess  the  feasibility  of  the  approach 
and  develop  the  required  techniques.  The 
considerations  include  the  materials  of  the  “tissue”  and 
fiducials,  3-D  imaging,  image  processing,  and  data 
reporting. 

2.1  Truth  Cube 

We  selected  a  2-part  silicone  rubber  (RTV6166, 
General  Electric  Co.)  because  of  its  proven  success  in 
simulating  tissues  [OttOO,  Wel99].  The  material  is  soft 
but  exhibits  linear  behavior  to  at  least  30%  strain.  To 
enable  tracking  of  the  internal  displacements  of  the 
cube  we  embedded  small  beads  as  fidueials  that  would 
readily  appear  in  CT  scans  but  not  significantly  alter 
the  material  properties  of  the  silicone.  The  choice  of 
bead  material  and  size  was  motivated  both  by  the  need 
for  high  contrast  without  creating  imaging  artifacts 
[Str98];  this  implies  a  specific  gravity  well  above 
silicone  (0.98)  but  much  less  than  steel  (~8).  To  enable 
accurate  segmentation  and  position  estimation,  the 
minimum  diameter  of  the  beads  had  to  be  slightly 
larger  than  the  distance  between  two  successive 
scanning  planes  (1.25  mm)  but  small  enough  to  avoid 
compromising  the  material  properties  of  the  silicone. 
Imaging  tests  revealed  that  Teflon  beads  (specific 
gravity  2.3)  with  a  diameter  of  L58  mm  fulfilled  all  of 
these  criteria. 

The  mold  for  the  cube  had  removable  sides,  with  the 
bottom  serving  as  a  rigid  permanent  mounting  plate  and 
was  roughened  and  dimpled  to  ensure  complete 
adhesion  of  the  silicone.  The  two-part  silicone  rubber 
was  mixed  in  a  30:70  proportion  to  obtain  material 
characteristics  similar  to  that  of  soft  tissue.  A  layer 
1  cm  deep  was  poured,  the  mold  was  placed  in  a 
vacuum  chamber  to  remove  air  bubbles,  and  then 
allowed  to  set  on  a  leveled  surface  (-3.5  hours)  before 
positioning  the  beads  in  a  7  x  7  matrix  spaced  1  cm 
apart  (Figure  1).  The  beads  were  laid  on  the  rubber 
using  a  matrix  of  bead  positioning  tubes  and  allowed  to 
set  into  the  rubber  for  at  least  another  hour  before  the 
next  layer  was  poured.  The  end  result  produced  an 
8cm3  silicone  rubber  cube  with  7  rows  and  7  columns 
of  Teflon®  beads  in  7  layers  spaced  1cm  apart 
(Figure  2). 

To  ensure  that  the  material  properties  of  the  silicone 
rubber  were  not  affected  by  the  addition  of  the  343 
beads,  we  performed  indentation  tests  on  a  sample  of 


Figure  1 :  (left)  Mold  and  bead  positioner  and  (right) 
resulting  “Truth  Cube” 


1.58  mm  diam 


Figure  2:  Dimensional  drawing  of  “Truth  Cube” 

the  silicone  rubber  alone,  on  a  sample  with  sparse  bead 
density,  and  on  a  sample  with  1  cm  bead  spacing. 
Results  indicated  that  the  stiffness  of  the  sample  with 
the  highest  bead  density  (780  N/m)  changed  2%,  less 
than  the  4.4%  standard  deviation  of  measurements  of 
the  plain  sample.  We  thus  assume  that  the  addition  of 
the  beads  had  a  negligible  effect  on  the  material 
properties  of  the  silicone. 

A  small  strain  (<4%)  uniaxial  compression  test  was 
performed  on  the  final  cube  to  obtain  material  property 
characteristics.  Results  indicate  a  Young’s  modulus  of 
25.8  kPa.  We  assume  that  Poisson’s  ratio  is  near  0.5  for 
this  polymer;  volumetric  date  from  the  loaded  state 
below  showed  minor  compression,  but  this  difference 
was  less  than  the  measurement  accuracy. 
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2.2  Experimental  Set-up 

The  test  set-up  for  loading  the  cube  must  apply 
carefully  controlled  boundary  conditions  during  the  CT 
imaging  process.  The  fixture  base  was  constructed  of 
25.4  mm  thick  PVC  (specific  gravity  1.37)  to  span  the 
width  of  the  scan  table  for  support  (Figure  3).  A  pocket 
that  tightly  fit  the  rubber  cube’s  base  plate  was  cut  into 
the  fixture  base.  To  ensure  accurate  positioning  of  the 
cube  each  time  it  was  tested,  various  reference  features 
were  incorporated  into  the  test  fixture.  The  fixture  base 
also  included  registration  markers  consisting  of 
machined  holes  and  grooves  in  “X”  patterns,  and  Delrin 
spheres  (specific  gravity  1.4)  positioned  at  the  cube 
base  comers  and  on  top  of  a  central  post  (visible  in 
Figure  4A,  B). 


Pivot  and 
vertical 


Figure  3:  Side  view  of  the  experimental  set-up.  L  is  the 
length  from  the  loading  points  to  the  center  of  the  cube, 
dz  is  the  vertical  motion,  and  m  is  the  mass  of  the  counter 
weights. 

The  initial  loading  test  reported  here,  uniaxial 
compression,  was  accomplished  with  a  25.4  mm  thick 
acrylic  plate  (specific  gravity  1.17),  This  was  attached 
by  a  low-friction  pivot  to  a  vertical  support.  The  plate 
was  loaded  by  metal  weights  applied  to  the  end 
opposite  the  pivot  (in  addition  to  the  weight  of  the  plate 
itself).  The  design  confined  all  metal  parts  (the  vertical 
support  rods  and  weights)  to  the  ends  of  the  apparatus, 
so  that  they  did  not  interfere  with  CT  scanning  of  the 
cube  in  the  middle  of  the  apparatus.  A  linear  dial 
indicator  measured  the  distance  the  plate  was 
compressed, 

2.3  Experimental  Protocol 

A  first  scan  of  the  experimental  setup  was  performed 
while  the  cube  was  unloaded.  The  flat  compression 
plate  was  then  held  level  as  it  was  lowered  onto  the 
oiled  top  surface  of  the  cube  to  a  set  displacement. 
With  the  pivot  end  of  the  compression  plate  held  in 
place,  masses  were  added  to  a  pocket  in  the  other  end 
until  equilibrium  was  established  as  determined  by 
level  indicators.  Three  loading  conditions  were 
scanned:  4  mm,  10  mm,  and  14,6  mm  displacements 


producing  5.0%,  12.5%,  and  18,25%  nominal  strain 
respectively. 

Once  loaded,  the  cube  was  scanned  in  a  General 
Electric  LightSpeed  CT  scanner.  Volumetric  images  of 
the  “truth  cube”  were  then  obtained  using  the  following 
scanner  settings:  120kV,  180mA,  standard 

reconstruction  type,  and  high  quality  scan  mode.  The 
field  of  view  was  of  190  mm,  producing  a  voxel  size  of 
0.37109  x  0.37109  x  1.25  mm3.  Each  CT  scan  was 
composed  of  99  slices  (512  x  512  pixels,  16  bits) 
exported  in  DICOM  format. 

2.4  Image  Processing 

Commercial  software  (3D-Doctor,  Able  Software 
Corp.,  Lexington,  MA)  was  used  for  manipulating  die 
volumetric  images.  First,  the  minor  artifacts  present  In 
the  lower  part  of  the  cube  (due  to  the  PVC  support 
plate)  were  filtered.  Second,  fiducial  beads  in  each  of 
the  scans  were  segmented  from  the  gel  with 
thresholding;  this  was  straightforward  due  to  the 
relatively  good  contrast  between  the  beads  and  the 
silicone  gel,  and  all  343  beads  were  readily  extracted. 
Surface  models  were  then  calculated  in  a  format 
suitable  for  further  processing. 

The  beads  in  the  silicone  cube  have  a  size  close  to  the 
size  of  a  voxel,  so  the  three-dimensional  segmentation 
process  only  provided  an  approximation  of  the  shape 
and  location  of  the  beads.  To  further  analyze  the  data  it 
was  necessary  to  accurately  determine  the  location  of 
the  center  of  each  bead.  We  implemented  a  least 
squares  algorithm  that  fit  a  sphere  to  the  set  of  points 
defined  by  the  vertices  of  the  segmented  beads.  This 
algorithm  is  part  of  visualization  software  that  visually 
assesses  the  quality  of  the  fit,  saves  the  information  to 
files,  and  displays  polygonal  models  and  vector  fields. 

Computation  of  the  displacement  vectors  of  the  beads 
between  various  loading  conditions  required  the 
solution  of  a  relatively  complex  matching  problem.  We 
therefore  used  a  manual  matching  approach  for  this 
first  set  of  experiments.  Finally,  the  outermost  three- 
dimensional  surfaces  of  the  cube  at  each  compression 
stage  were  extracted  from  the  images.  Registration 
marks  were  also  segmented  in  order  to  evaluate 
potential  motion  of  the  experimental  setup  between 
successive  scans, 

2.5  Finite  Element  Model  Comparison 

A  finite  element  model  of  the  truth  cube  using 
commercial  FEM  software  (ABAQUS  v.  5.8,  Hibbitt, 
Karlsson  Sc  Sorensen,  Inc.,  Pawtucket,  Rhode  Island) 
was  computed  to  compare  with  the  experimental 
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results.  A  mesh  was  manually  generated  consisting  of 
10  mm,  8-node,  solid  cube  linear  elements.  The  cube’s 
material  properties  (as  measured  in  [OttOl])  are 
isotropic  and  linear  with  a  Young’s  modulus  of  15.3 
kPa  and  an  assumed  Poisson’s  ratio  of  0,499.  Large 
deformation  (i.e.  nonlinear)  analysis  was  employed 
using  the  boundary  conditions  and  geometry  extracted 
from  the  CT  scans  for  the  unloaded  and  18,25% 
nominal  strain  cases, 

3.  Results 

Figure  4  illustrates  the  processing  of  the  data, 
beginning  with  segmented  CT  data  for  the  x-z  slice 
through  the  center  of  the  cube,  in  both  unloaded  (Figure 
4A)  and  18.25%  strain  (Figure  4B)  cases.  Figure  4C 
shows  the  bead  displacements  throughout  the  cube 
combined  for  all  four  stain  levels,  while  Figure  4D 
represents  the  displacement  field  for  the  central  x-z 
slice.  Reconstructions  of  the  outer  surface  of  the  cube 
for  each  strain  state  is  shown  in  Figure  4E,  The 
estimated  error  of  the  bead  locations  is  2  mm. 

The  results  from  the  FEM  model  are  shown  in 
Figure  5  A,  in  the  form  of  a  plot  for  the  x-z  slice  through 


the  center  of  the  cube  with  displacement  magnitude 
coded  as  color.  Figure  5B  shows  the  FEM  and 
experimental  displacement  vectors  for  each  bead  in  this 
slice.  The  coordinates  of  a  sample  of  10  of  the  beads 
are  presented  in  Table  1,  together  with  their 
displacements  from  the  FEM  and  experimental 
measurements,  and  the  differences  between  these 
estimates.  Full  results  are  available  on  the  project  web 
site,  http://biorobotics.harvard.edu/truthcube/. 

4,  Discussion 

The  results  presented  here  demonstrate  the  potential  of 
this  technique  for  establishing  physical  standards  for 
validating  real  time  soft  tissue  simulations.  The 
construction  of  the  cube  was  straightforward  and 
resulted  in  a  structure  with  a  modulus  in  the  range  of 
soft  organs.  Imaging  under  strains  of  over  18%  was 
readily  accomplished.  The  external  surface  and  the 
fiducial  beads  were  easily  segmented,  which  permitted 
calculation  of  the  displacement  fields  throughout  the 
volume.  The  result  is  a  displacement  data  set  with  well- 
characterized  material  properties  and  boundary 
conditions. 


Figure  4:  The  CT  image  of  the  eentral  x-z  plane  of  the  "Truth  Cube”  after  segmentation  in  both  the  (A)  unloaded  and  (B)  maximally 
compressed  states;  the  upper  registration  sphere  is  noticeable  in  these  images,  (C)  3_d  perspective  view  of  the  centroids  of  best-fit  spheres 
for  all  343  beads  in  all  strain  states  (blue  =  unloaded,  green  =  5%,  red  =  12.5  %,  yellow  =  18.25%  strain).  (D)  Trajectory  of  the  beads  in 
the  x-z  central  plane  of  (C).  (E)  Outer  surfaces  of  the  cube  in  its  4  strain  states  (purple  -  unloaded,  red  =  18,25%). 
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Figure  5:  A  comparison  of  the  finite  element  displacement  field  and  the  CT  experimental  data  displacement  field  for  the  central 
XZ  plane.  (A)  A  picture  of  the  displacement  of  the  uncompressed  (black)  and  compressed  (colored)  mesh.  (B)  A  comparison  of 
the  finite  element  results  (thin  black  lines)  and  experimental  results  (grey-scale  lines). 
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3.1 

0.6 

14.9 

3.0 

-9.0 

3.0 

13.0 

-0,6 

0.0 

0.6 

-0.6 

0.4 

1.2 

0.7 

1.0 

3.0 

13.0 

0.0 

0.0 

0,7 

0,1 

0.2 

1,1 

0.4 

11,0 

3.0 

13.0 

0.6 

0.0 

0,6 

0.7 

0.2 

1,3 

0.7 

21,0 

3,0 

13.0 

1.3 

0.0 

0,8 

1.8 

0.1 

1,2 

0.6 

31.0 

3.0 

13.0 

2,1 

0.0 

0,6 

3.0 

0.0 

1.8 

1,5 

Table  1:  A  representative  set  of  10  bead  locations  at  rest  and  after  18.25%  strain  (14.6  mm  compression)  evaluated  by  both 
FEM  and  experimental  data.  The  final  column  is  the  difference  in  bead  location  between  the  two  evaluations. 


These  results  are  readily  compared  with  the  output  of 
computational  models.  The  simple  FEM  calculation 
presented  here  used  a  minimal  number  of  elements  to 
reflect  some  of  the  speed-vs-accuraey  tradeoffs  of  real¬ 
time  tissue  simulations.  The  maximum  differences  with 
the  experimental  results  were  on  the  order  of  3.5  mm 
for  a  displacement  of  about  15  mm.  This  is  a  relatively 
small  difference,  as  might  be  expected  for  these  simple 
geometric  and  loading  conditions,  and  considering  the 
estimated  2  mm  accuracy  of  the  measurement.  Some 
portion  of  the  difference  may  also  be  due  to  inexact 
specification  of  the  boundary  conditions  on  the  top  of 
the  cube  during  compression.  Nonetheless,  refinement 
of  the  model  (e.g.  smaller  mesh  size)  should  bring 
closer  agreement  to  the  experiment.  For  more  complex 
geometries,  materials,  and  loading  conditions,  much 
larger  differences  can  be  expected  between  FEM  results 
and  reality,  which  this  type  of  experimental  data  can 
help  quantify. 


Several  refinements  to  the  measurement  techniques  are 
indicated  by  these  results.  First,  we  have  only  estimated 
the  accuracy  of  the  bead  positions  in  this  initial  dataset. 
More  precise  accuracy  specifications  can  be  obtained 
through  an  analysis  of  the  imaging  and  image 
processing  procedures,  and  through  reproducibility 
measurements.  The  raw  scan  data  will  be  available  on 
the  web  site,  so  other  users  are  also  welcome  to  develop 
improved  image  processing  algorithms  for  this  purpose. 
Another  area  for  better  characterization  is  the  top 
surface  boundary  condition,  where  oiling  the  cube 
served  to  limit  and  regularize  friction,  but  where  a 
“frictionless”  condition  clearly  did  not  apply. 
Fortunately,  the  contact  area  is  readily  measured  during 
the  experiment  and  in  the  processed  CT  images,  so 
additional  boundary  information  is  available. 

4,1  Future  Work 

Our  immediate  plans  are  to  complete  tests  of  the  truth 
cube,  with  a  goal  of  improving  the  accuracy  of  the 
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fiducial  position  estimation  process  and  the  boundary 
conditions.  We  will  also  develop  additional  loading 
capabilities,  including  indentation,  shear,  and  torsion, 
to  better  reflect  the  range  of  surgical  manipulations. 

While  this  initial  study  demonstrated  the  feasibility  of 
experimental  measurement  of  volumetric  displacement 
fields  for  soft  materials  under  large  strains,  the  truth 
cube  itself  is  of  limited  value  for  validating  soft  tissue 
models.  The  cube  has  regular  geometry  and  material 
properties  that  are  helpful  for  the  development  process, 
but  fast  tissue  simulations  must  deal  with  conditions 
that  are  vastly  different,  involving  very  large 
deformations,  irregular  shapes  and  complex  materials. 
These  are  also  the  conditions  where  validation  using 
traditional  FEM  techniques  is  most  problematic. 

The  next  phase  of  this  project  is  thus  the  measurement 
of  biological  tissue  volumetric  displacements,  in  the 
form  of  an  entire  organ.  Preliminary  plans  call  for 
measurement  of  a  porcine  liver  immediately  post 
mortem.  Fiducials  will  be  inserted  via  needles 
throughout  the  parenchyma,  followed  by  the  imaging 
and  image  processing  procedures  reported  here.  In 
addition,  a  cast  of  the  external  shape  of  the  liver  will 
enable  the  creation  of  a  silicone  rubber  model  with 
bead  fiducials  where  the  material  properties  will  be 
well-characterized.  Simulation  developers  can  then 
separately  consider  issues  of  large  deformation  and 
material  complexities. 

While  these  biological  material  tests  have  obvious 
advantages,  the  material  properties  will  not  be  well- 
characterized  as  in  the  rubber  tests.  The  identification 
of  these  material  properties  is  an  active  area  of 
research,  and  the  data  sets  may  be  useful  for  solving  the 
inverse  problem  of  estimating  these  material  properties 
from  the  displacement  data.  In  any  case,  the  data  will 
provide  a  source  of  information  for  the  development  of 
fast  simulation  systems. 

As  we  proceed,  we  welcome  comments  and 
suggestions  on  all  aspects  of  the  project,  from  tissue 
selection  and  loading  conditions  to  experimental 
methods  and  data  analysis.  Further  information  is 
provided  on  the  web  site  at: 

http://biorobotics.harvard.edu/truthcube.  Contributions 
of  a  broad  cross-section  in  the  tissue  modeling 
community  will  help  ensure  the  development  of  a 
robust  physical  standard  for  the  benchmarking  of  real¬ 
time  soft  tissue  simulations, 
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